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Abstract —An important enabler towards the successful deploy¬ 
ment of any new element/feature to the cellular network is the 
investigation and characterization of the operational conditions 
where its introduction will enhance performance. Even though 
there has been significant research activity on the potential of 
devlce-to-devlce (D2D) communications, there are currently no 
clear indications of whether D2D communications are actually 
able to provide benefits for a wide range of operational conditions, 
thus justifying their Introduction to the system. This paper 
attempts to fill this gap by taking a stochastic geometry approach 
on characterizing the set (region) of operational conditions for 
which D2D communications enhance performance in terms of 
average user rate. For the practically Interesting case of a 
heavy loaded network, the operational region is provided in 
closed form as a function of a variety of parameters such as 
maximum D2D link distances and user densities, reflecting a wide 
range of operational conditions (points). It is shown that under 
the appropriate deployment scheme, D2D communications can 
indeed be beneficial not only for the usually considered regime of 
“proximal communications” but to a wide range of operational 
conditions that include D2D link distances comparable to the 
distance to the cellular access point and considerably large user 
densities. 

Index Terms —cellular network, channel-access-mode selection, 
D2D communications, optimal system design, stochastic geometry 


I. Introduction 

The exponential growth of traffic that the cellular networks 
are currently experiencing has driven academia and industry 
to consider new channel access methods in addition to the 
conventional, infrastructure-based access via cellular nodes 
(access points). One of these methods that have recently 
attracted a lot of attention is device-to-device (D2D) communi¬ 
cations as a means for offloading cellular traffic as well as op¬ 
portunistically exploiting link proximity ||2|. Even though the 
potential of D2D communications cannot be denied, even more 
so with the ever increasing number and computational/storage 
capabilities of devices, their introduction as part of the cellular 
network introduces new challenges El with respect to (w.r.t.) 
resource sharing between cellular and D2D access modes. 

Specihcation of (optimal) resource sharing is of great im¬ 
portance, nevertheless, it is not alone sufficient to justify 
the introduction/adoption of a D2D mode as an “add-on” to 
the (conventional) cellular network, even if the corresponding 
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analysis shows great potential under certain operational con¬ 
ditions. For example, it is obvious that exploitation of D2D 
links with “small” distance compared to the distance from the 
closest access point is benehcial, which actually provides the 
incentive to investigate D2D communications in the hrst place. 
However, the vague notion of “closeness” (proximity) must be 
converted to actual figures in order to assess the merits of D2D 
communications, since it is reasonable to wonder whether the 
proximity requirements of D2D communications are actually 
so strict that render the benefits of D2D communications non¬ 
existent for realistic/practical conditions. Similar concerns can 
be posed w.r.t. other relevant parameters such as (relative) 
densities of concurrent D2D and cellular links. 

It becomes clear that a detailed investigation and char¬ 
acterization of the range of operational conditions where 
introduction of a D2D mode enhances cellular performance 
is a procedure of critical importance. These operational con¬ 
ditions are typically described by quantities that the system 
designer/operator has no control over, such as density of users 
and channel conditions, as well as quantities chosen at a 
pre-deployment system design stage such as maximum link 
distance for which D2D transmission can (potentially) be es¬ 
tablished. Characterization of the D2D operational region will 
determine whether the introduction of D2D communications 
is actually beneficial for a (hopefully) large dynamic range 
of use/business cases and, therefore, plays a key role in the 
adoption of a D2D mode as part of the cellular network. 

A. Related Work 

There has been considerable work on resource sharing 
for D2D communications under an optimization modeling 
framework (see Q for a recent literature review). Even 
though providing optimal solutions, this framework is bound to 
consider finite-area networks, and, therefore, provides limited 
insight for large-scale network deployments. Toward this end, 
a second line of works has recently emerged, attempting to 
describe the merits of D2D communications based on tools 
from stochastic geometry 0, which have been successfully 
applied for the analysis of ad-hoc and cellular networks lb], 

Q. 

Out of the relatively small number of notable publications 
considering analysis and design of D2D communications un¬ 
der a stochastic geometry framework, the most relevant with 
this work are 0-111. An overlay D2D network deployment, 
i.e., where cellular and D2D transmissions are performed in 
orthogonal (non-overlapping) bands of the available spectrum, 
is considered in 0. The distribution of transmit powers is 
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investigated under a power control scheme, however, without 
any considerations regarding spectrum sharing between cellu¬ 
lar and D2D networks. An underlay D2D network deployment, 
i.e., where cellular and D2D transmissions are performed in 
the same bandwidth, is considered in i). A simple on-off 
power control scheme for D2D transmissions is proposed in a 
setup considering only a single cell serving one uplink user. 
More comprehensive studies on the co-existence of large-scale 
deployed cellular and D2D networks are pursued in lfT0l - lfT4l . 
In noi-iiiii, underlay or overlay D2D communications are 
considered along with an uplink or downlink cellular network. 
System operation aspects such as resource sharing among 
networks and mode selection according to various criteria 
are investigated. In ini, system performance in terms of 
spectral efficiency (average rate per channel use) is considered 
under a scheme providing exclusion regions to protect cellular 
users from D2D-generated interference, whereas in llT4l . the 
transmission capacity, i.e., density of cellular and D2D links 
with signal-to-interference ratio above a certain threshold, is 
considered. All the above works, even though they provide 
valuable insights on the potential of introducing a D2D mode 
to the cellular system, do so by considering only a limited 
number of (artificially selected) operational points (scenarios) 
without any clear indication of the region of operational 
points where D2D communications are beneficial to system 
performance. 

B. Contributions 

This paper considers D2D communications deployed as 
overlay or underlay to the downlink-dedicated cellular band¬ 
width. The main goal is to provide a mathematically rigorous 
specification of the operational region where introduction of a 
D2D mode enhances cellular network performance in terms of 
average user rate. The operational region is described in terms 
of the following parameters (stated in no particular order) 
whose possible combinations of values can be mapped to a 
wide variety of use/business cases of potential interest: 

• density of users, 

• density of cellular access points, 

• propagation conditions, 

> signal-to-interference-ratio (SIR) threshold for successful 
communications, and 

> maximum D2D link distance. 

A number of D2D deployment schemes that trade-off 
complexity for performance are investigated, characterized 
by whether (non-trivial) mode selection and/or D2D channel 
access procedures are employed. For the mode selection 
procedure, two alternatives are considered depending on the 
availability of D2D link distance information. Based on ac¬ 
curate analytical expressions for the average rates achieved 
by the cellular and D2D links, a tractable expression for 
the average user rate of a D2D-enabled cellular network is 
obtained. For the practically interesting case of a heavy loaded 
system (large user density), the operational regions for each 
D2D deployment scheme are obtained in closed form, allowing 
for investigating the effect on system performance of the two 
most important parameters related to D2D use/business cases. 


TABLE I 

Summary of Key Notation 


Notation 

Description 

R, R+, 0 

set of reals, set positive reals, empty set 

P(.),E(.),I(-) 

probability measure, expectation operator, indicator 
function 

Isl 

Euclidean norm of x G 

•^0,9 ^di 

densities of APs, D-UEs, and C-UEs, respectively 

rd 

distance of a random (potential) D2D link (r^ £ 

(0; ^(i,max]) 

Td.th 

threshold distance used by the distance-based mode 
selection procedure (rd^h G (0, rd^max]) 

Pa 

(common) AP transmit power 

P 

probability of D-UE selecting D2D mode {p £ (0,1]) 

<? 

probability of transmission of a D2D link (q £ (0,1]) 

a 

path loss exponent (a > 2) 

K, Ko 

number of cellular UEs within a random cell and the 
cell containing the origin, respectively 

Vc 

portion of bandwidth dedicated to cellular transmissions 
(Pc e (0,1]) 

K 

p{e) 

^ (27r/a)/sin(27r/a) 

- fe-2/. 

7 

= 1, 2 for probabilistic and distance-based mode selec¬ 
tion, respectively 

^0 

SIR threshold below which communication is consid¬ 
ered unsuccessful (0o > 0) 

Fop, Vi 

system operational point, system design parameters 

R, Rc, Rd 

rate of typical link, typical cellular link, and typical D2D 
link, respectively 

RnoDlD 

rate of typical cellular link when the network does not 
support D2D transmissions 

7?^ 

Operational region of D2D communications using 
Scheme i £ {1, 2, 3, 4} 


namely, maximum D2D link distance and density of users 
(devices) capable for D2D communications. 

Important outcomes of this analysis include the following: 

1) D2D communications can provide benefits for opera¬ 
tional conditions of practical interest including D2D link 
distances of the order of the distance to the closest access 
point. 

2) Employing both a mode selection and a D2D channel 
access procedure does not necessarily provide significant 
gains compared to employing one of them only. 

3) Exploitation of D2D link distance information alone for 
mode selection purposes is sufficient for D2D commu¬ 
nications to be beneficial for any operational point. 

II. System Model 

A. System Nodes Description and Channel Access Mecha¬ 
nisms 

A D2D-enabled cellular system is considered that allows for 
overlay or underlay downlink-inband D2D transmissions, i.e., 
transmissions are performed either in cellular or D2D mode. 
In particular, the system serves two types of user equipments 
(UEs) E). The first type, referred to in the following as C- 
UEs, request data from sources that do not have the ability 
to establish a direct communication link with the requesting 
C-UE, due to being located far away and/or not having 
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transmit capabilities. C-UEs can only utilize downlink cellular 
transmissions in order to obtain their data. The second type 
of UEs, referred to in the following as D-UEs, request data 
from sources that are physically proximal and have transmit 
capabilities, providing D-UEs the flexibility to establish either 
a D2D or a dual-hop (uplink and downlink) cellular link for 
communication. Eull buffers are assumed for all communica¬ 
tion links. 

The positions of C-UEs and D-UEs are modeled as inde¬ 
pendent homogeneous Poisson point processes (HPPPs) <i)c, 
C of densities Ac, Xd > 0, respectively, whereas 
the position of the source for a D-UE located at cc G <i>d 
is modeled as an independent random variable uniformly 
distributed within the closed ball {y S : |t/ —x| < rd^max}- 
By application of the displacement theorem for HPPPs 0, it 
can be shown that the sources of D-UEs also constitute an 
HPPP of density Xd which, with a slight abuse of notation, 
will be also denoted as The maximum possible D2D link 
distance rd,max > 0 provides an exact description of the notion 
of “proximity” and is determined by the relevant D2D use 
cases and/or imposed by considerations regarding, e.g., power 
consumption or system performance. Note that parameters Ac, 
Ad and Vd^max are correlated, e.g., increasing rd,max can poten¬ 
tially “transform” some of the C-UEs into D-UEs. However, 
there are currently no well-established models available in 
the literature for modeling this correlation. Therefore, in this 
paper, these parameters are treated as independent variables, 
with the corresponding analysis easily extended under any 
model describing their interrelation. 

The cellular system infrastructure consists of access points 
(APs) randomly deployed according to an independent HPPP 
C R^ of density Aa > 0 and each UE employing cellular 
transmissions is served by its closest AP. In order to avoid 
intracell interference, time is (universally) divided into discrete 
slots of equal duration and UEs within a cell employing 
cellular transmissions are served by a time-division-multiple- 
access (TDMA) scheme with no priorities among them (round- 
robin scheduling). All active APs, i.e., APs with at least 
one served UE, transmit with a fixed power Pa, without any 
coordination/cooperation among them. 

In contrast to cellular communications where intracell in¬ 
terference is eliminated by centralized scheduling, a proba¬ 
bilistic channel access scheme fb), ifTSl is employed for D2D 
communications. Specifically, D2D communications utilize the 
same time slots as the cellular system (perfect synchronization 
is assumed between D2D and cellular slots in the underlay 
case), with the source of any established D2D link attempting 
to transmit in each slot with a probability q G (0,1]. Active 
D2D links employ open loop power control so that the effect 
of path loss, averaged over small-scale fading, is eliminated 
(large-scale path loss inversion) IfT^ . IfT^ . 

B. Mode Selection and D2D Network Deployment Schemes 

Since D-UEs have the flexibility to choose between two 
access methods, a mode selection procedure should be em¬ 
ployed. In this paper, the following two simple procedures are 
considered: 


Definition 1: Probabilistic mode selection: Each D-UE se¬ 
lects D2D mode by independently tossing a biased coin with 
bias pprob e ( 0 , 1 ]. 

Definition 2: Distance-based mode selection: Each D-UE 
selects D2D mode if and only if the distance Vd G (0, r^^max] 
from its source is smaller than a pre-defined threshold th G 

( 0 i ^d,inax]- 

Note that both mode selection procedures can be viewed 
as thinning operations 0 on with retention probability 
equal to pprob and P(rd < rd,th) = (rd,th/rd,max)^ for the 
probabilistic and distance-based case, respectively. This allows 
for a unified treatment of both mode selection schemes under 
a common notation p G (0,1] for their corresponding D2D 
mode selection probabilities. Parameters p and q constitute the 
so called D2D mode parameters and their (optimal) values will 
be shown to require knowledge of parameters such as AP and 
UE densities. Therefore, Ppiob, rd,th, and q are provided to the 
D-UEs via a broadcast mechanism. 

Since incorporation of mode selection and/or channel access 
procedures comes with implementation cost, e.g., for time 
slot synchronization or acquisition of D2D link distance in¬ 
formation, it is of interest to examine and compare the merits 
of each of the following D2D network deployment schemes, 
trading off flexibility on design of D2D mode parameters to 
complexity: 

• Scheme 1 (Baseline): p = 1, q = 1, 

• Scheme 2 (D-UEs employ D2D mode by default): p = 1, 
q G (0,1], 

• Scheme 3 (D2D links always active): p G (0,1], 9 = 1, 

• Scheme 4 (Most general case): (p,q) G (0,1]^. 

Scheme 3 has two versions, according to whether prob¬ 
abilistic or distance-based mode selection is employed, that 
will be denoted, when necessary, as 3-p and 3-d, respectively, 
with similar notation also used for Scheme 4. It is noted that 
Schemes 3-d and 4-p were first examined in Eol and ca, 
respectively, under different system model and/or performance 
metrics than this paper. Clearly, Scheme 4 is optimal under any 
performance criterion as it incorporates the other schemes as 
special cases. However, the implementation cost of operating 
under Scheme 4 with both p < 1 and q < 1 must be justified 
by a significant performance gain compared to the simpler 
schemes, an issue that will be examined in the following 
sections. 

Remark: Unless stated otherwise, results and definitions in 
the following sections correspond to operation under Scheme 
4. Translation to the other schemes can be made by trivially 
setting p and/or q equal to 1 . 

C. Signal-to-Interference Ratio 

The standard approach of conditioning on the existence of a 
(typical) UE located without loss of generality (w.l.o.g.) at the 
origin will be employed for the following analysis. Note that, 
by the properties of HPPP, this conditioning does not have 
any effect on the distribution of the other UEs 0. Treating 
interference as noise and neglecting the effect of thermal noise 
due to the system operating in the interference limited region, 
the SIR experienced by the typical UE when served by a 
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transmitter (either an AP or a proximal source) located at 
yo G equals 

{PaKVO e $a) + |2/ori(yo e $d)) ffyol2/or“ 

yeia.\vo v&^d\vo 

where C $„, C <i)(j denote the sets of positions of APs 
and D-UE sources, respectively, that access the channel at the 
considered time slot and bandwidth, gy is the channel gain 
corresponding to a transmitter located at y G tha U Xy is 
the position of a D-UE receiving data from its source located 
at y G and a > 2 is the path loss exponent. The channel 
gains are assumed independent and exponentially distributed 
(Rayleigh fading) with mean one. 

Note that w.l.o.g. the useful received power of a D-UE in 
D2D mode, averaged over channel fading, is normalized to 
one, which allows to treat Pa as the only variable that controls 
the relative interference levels of cellular and underlay D2D 
transmissions. Under this power control model, the transmit 
power of the source of a D2D link may become unrealistically 
large when D2D links of excessively large distances are 
established. In order to study the fundamental gains/limitations 
of employing a D2D mode as part of a cellular network w.r.t. 
to average user rate, no constraint on D2D maximum transmit 
power is set in this paper. However, as it will be shown 
later for operational scenarios of practical interest, significant 
performance gains are only provided by establishing D2D 
links of distances up to about the order of the distance of a D- 
UE from its closest AP, corresponding to power requirements 
similar to uplink cellular transmissions. 

III. Average Cellular and D2D Rates 

The SIR is a critical measure of performance as in¬ 
creased SIR translates to improved spectral efficiency (in 
bits/Hz/channel use). In particular, the commonly employed 
mapping SIR log 2 (l + 0o)I(SIR > Oq) will be used for 
determining the spectral efficiency of a link 0 , ini, m, 
where 0o > 0 is the SIR threshold below which communica¬ 
tions are considered unsuccessful. However, the actual rate 
achieved by a UE is also affected by the multiple access 
procedures, namely, TDMA for cellular links and probabilistic 
channel access for D2D links0 This section provides closed 
form expressions for the average rate of cellular and D2D links 
of the typical UE, which, in turn, will determine the more 
important metric of average UE rate that will be considered 
in the next section. 

A. SIR Distribution of Cellular and D2D Links 

The SIR distribution of cellular and D2D links, required 
to characterize the statistical properties of the corresponding 
spectral efficiencies, are provided in Lemmas 1-2. 

Lemma 1: The SIR distribution of the typical UE when 
employing a cellular link in a system with underlay D2D 

*The effect of issues such as signaling, e.g., for establishment of transmis¬ 
sion links, is ignored. 


communications is 


p(siR > e) 


1 


i + ¥{K>Q)p{e) 


qp'fKX^rj^ 

2Xa 


(A) 


2/a ’ 


( 2 ) 

for 6 > 0, where K is the number of cellular receivers (RXs), 
composed from C-UEs and D-UEs selecting cellular mode, 
within a random cell of the network, 7 = 1 , 2 for probabilistic 
and distance-based mode selection, respectively, and k, p{6) 
as defined in Table I. In case of overlay D2D communications, 
Eq. (| 2 ]l holds with the last term of the denominator removed. 

Proof: See Appendix A. ■ 

Lemma 2: The SIR distribution of the typical UE when 
employing a D2D link in a system with underlay D2D 
communications is 


P(SIR > 6») Ri exp[-At7r6i2/“((l/2)gp^Adr^ 

+ AaP(iT>0)PE“)], (3) 

for 6 > 0, where 7 = 1, 2 for probabilistic and distance- 
based mode selection, respectively. In case of overlay D2D 
communications, (O holds with the last term inside parenthesis 
of the exponential removed. 

Proof: Proof follows the lines of ca and is omitted. ■ 
Remark: The expressions of Lemmas 1 and 2 are approx¬ 
imate since they are based on the assumption that is 
generated by a thinning of with a retention probability 
P(iT > 0) for each x G independent of everything else, 
which was shown in El to be a very good approximation. 
Eor the case when ¥{K > 0) = 1, i.e., when all APs serve at 
least one UE each, and, therefore, = $0 almost surely, © 
and © are exact ifTTl . Note that an approximate closed form 
expression for F{K > 0) with very good accuracy is available 
ifl^ . namely, 

P(iT > 0) 1 - (^1 -f . (4) 


B. Average Rate of Cellular Link 

Exact computation of the rate achieved by the typical UE 
when employing cellular communications is difficult since 
there exist cases when the end-to-end link will consist not 
only of a downlink, but also of an uplink hop, e.g., when 
the typical UE is a D-UE selecting cellular mode. In order to 
simplify the problem, the following assumption is made. 

Assumption 1: The cellular system is downlink limited in 
terms of user rate. 

Assumption 1 is reasonable when the number of cellular 
RXs requesting data from the network core (e.g., an internet 
server) is (much) larger than the number of cellular RXs 
requesting data from another UE, with the end-to-end cellular 
rate essentially limited by the downlink TDMA process. 
Towards describing the effect of TDMA on the average rate of 
the typical UE employing cellular transmissions, the following 
lemma will be of use. 

Lemma 3: Let Kq and K denote the number of cellular 
RXs positioned within the cell containing the origin and any 
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Other randomly selected cell, respectively. It holds 


/ 1 \ ^ AaP(K > 0 ) 

\ Ko + 1 / Ac + (1 — p) Ad 


(5) 


Proof: See Appendix B. ■ 

Note that the term 1/{Kq + 1) < 1 of the above Lemma 
represents the rate reduction of the typical cellular UE rate 
due to resource sharing via TDMA. 

Remark: Previous works considering the effect of TDMA 
either relied on numerical evaluations based on complicated, 
approximate expressions for the distribution of Kq m, eqi, 
or employed heuristic approximations of E(l/(Aro + 1 )) with 
the corresponding expressions equal to Eq. © with V{K > 0) 
replaced by 1 El or 7/9 El, ESI- The exact expression of 
Lemma 3 shows that these approximations do not hold in the 
general case. 

The average rate of the typical UE when employing cellular 
transmissions, normalized by the total bandwidth dedicated for 
downlink cellular and D2D transmissions can now be obtained 
as follows. 

Proposition 1: The average normalized rate achieved by the 
typical UE when employing cellular mode equals 


Rc = PcE 


1 


7^0 + 1 
77cAcP(iT > 0 ) 


Ac + (1 ~ p)Ad 


I(SIR> 0 o)jlog 2 (l + 0 o) ( 6 ) 

P(SIR>0o)log2(l + 0o), (7) 


in bits/s/Hz, where rjc is the portion of the total bandwidth 
where downlink cellular transmissions take place (0 < ? 7 c < 1 
and rjc = 1 for overlay and underlay D2D communications, 
respectively) and P(SIR > Oq), ¥(K > 0) are given by Eqs. 
I© and I©, respectively. 

Proof: By ignoring the correlation between Kq and 
SIR due to their common dependence on can be 

approximated as 


Rc 


ric^ P(SIR > Oq) log 2 (l + Oq). ( 8 ) 


The result then follows by application of Lemmas 1 and 3. ■ 
Remark: The expression on the right hand side of © was 
used as the actual definition of cellular rate in previous works, 
e.g., El, ED, implicitly suggesting that Kq and SIR are 
uncorrelated (or even independent) which does not hold. As 
will be shown in Sec. VII, this approximation is very accurate, 
especially for the overlay D2D case. 


C. Average Rate of D2D Link 

Considering the same SIR threshold 0 q as for cellular 
communications, the average normalized D2D link rate, Rd, 
is as follows (compare with the corresponding expression of 
Rc). 

Lemma 4: The average normalized rate achieved by the 
typical UE when employing overlay D2D mode equals 

Rd = (1 —Pc)E [I(D2D link active)I(SIR > 0o)] log 2 (l + ^o) 
= (l-pe)gP(SIR > Oq) log 2 (l + Oq) (9) 

in bits/s/Hz, with rjc as in Proposition 1 and P(SIR > 9 q) 
as given in Eq. ©. The rate expression for underlay D2D 


transmissions is the same as overlay with the term 1 — rjc 
replaced by 1 . 


IV. Average UE Rate AND System Design 
Considerations 


In this paper, incorporation of D2D communications is 
considered as a means for enhancing the performance of the 
conventional cellular system. Towards this end, the average 
normalized rate achieved by the typical UE without a-priori 
information of its type (C-UE or D-UE) is considered as the 
system metric of interest, equal to M, ED 

R = [pRci + (1 - p)Rc) (10) 

■Ac “r Ac “r Ad 


in bits/s/Hz, with Rc, Rd, the average rate of the cellular and 
D2D links, respectively, as described in Sec. III. Note that the 
rate achieved by the typical UE when a D2D mode option is 
not available, denoted as i?noD 2 D, is a degenerate case of (O 
for p = 0 and rye = 1 , and equals 


7 ?noD 2 D 


AaP(i^ >O)log2(l + 0o) 

(A, + Ad)(l+P(iT>O)p(0o))’ 


( 11 ) 


as can be easily verified using the results of Sec. III. 

Average rate R depends on a large number of parameters 
(through Rc and Rd) which are not shown explicitly for ease 
of notation. These parameters can be partitioned into two 
categories. The first category describes the system operational 
point and is represented by the tuple 


Rop — (Ac, Ad, Aa, 


K, 00 , Ct) , 


which consists of parameters that the system operator has no 
control over, namely, Ac, Ad, and a, and parameters which 
are chosen at system design stage and remain fixed during 
operation, namely, Aa, rd,max, and 0 q. The second category of 
parameters is represented by the tuple 

T’d = {p,q,r]c,Pa), 


which contains the system design parameters that can be 
selected by the operator in accordance to Pop- Note that 
specification of Pa is irrelevant for system design with overlay 
D2D, whereas 77 c = 1 by default with underlay D2D. There¬ 
fore, Pd essentially consists of the two D2D mode parameters 
and one parameter reflecting the resource sharing between 
cellular and overlay/underlay D2D networks. 

Since the D2D mode is considered as an enhancement 
option to the conventional cellular network, it is natural to 
impose that its introduction should not degrade the perfor¬ 
mance experienced by UEs employing cellular transmissions 
0 , i). This can be reflected by the following restriction on 
the values of Pd. 

Constraint 1: Pd is allowed to take values for which the 
average normalized rate of the cellular links is unaffected by 
the introduction of a D2D mode, i.e., Rc — i?noD 2 D- 

Constraint 1 is convenient as it allows for eliminating one 
element of Pd, thus reducing the design space and simplifying 
the problem. In particular, by employing Eqs. I© and ET) . 
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the resource sharing design parameters Pa and r]c can be 
determined in closed form as functions of p and q, namely. 


/ {^c + {l-p)Xd)Krl^^^qp^ 

2Xa{l + V{K >0)p{9o)) J 


for the underlay D2D case, with 7=1,2 for probabilistic and 
distance-based mode selection, respectively, and 


ric = l- 


P>^d 

Ac + Ad ’ 


(13) 


for the overlay D2D case, irrespective of the mode selection 
scheme. Notice how the complicated interference environment 
that a UE on cellular mode experiences in the underlay D2D 
case is reflected on the expression for Pa, which increases 
as in order to compensate for the correspondingly 

increasing transmit power of D2D links. In contrast, the 
expression for rjc in the overlay case is much simpler, having 
the intuitive interpretation that bandwidth is proportionally 
partitioned according to the relative densities of cellular and 
D2D links. 

Specification of the resource sharing parameters by Eqs. 
(fT2li and (fTJt leaves the D2D mode parameters as the only 
variable elements of Pd- Their optimal values 


(p*,q*)=a.Tg max R (14) 

(P,9)6(0,1]2 

can then in principle be found by a numerical search using 
the analytical formulas of the previous section for any oper¬ 
ational point T^op of interest. Another quantity of significant 
importance is the D2D operational region defined as follows. 


Definition 3: The D2D operational region, P-mt), is the set 
of system operational points for which incorporation of a 
D2D mode can increase the average UE rate provided by the 
(conventional) cellular network, i.e., 

"^020 = {Pop ■ 3 (p, q) G (0, 1]^ for which R{Pop) > i?noD2D}- 

(15) 

Note that 7 ?.d2d includes operational points where an arbitrar¬ 
ily small rate increase is achieved. These points, even though 
of small interest in that sense, may still be considered as 
beneficial due to implicit gains achieved w.r.t. other metrics, 
e.g., signaling levels, latency. 


V. Average UE Rate and Ordering of D2D 
Deployment Schemes for Heavy Loaded System 

A. Average UE Rate for Heavy Loaded System 

The closed form expressions of the previous sections allow 
for efficiently computing {p*,q*) for any P^^ of interest. On 
the other hand, obtaining 7 ?.d2d numerically is a rather compli¬ 
cated task as the domains of (p, q) and P^^ are uncountable. 
In order to obtain closed form expressions for (p*, q* ) and 
Pd 2 d the following assumption will be employed. 

Assumption 2: The density Ac of C-UEs is sufficiently large 
so that F{K > 0) Ri 1 (irrespective of the density Xd of D- 
UEs), where K is the number of cellular RXs in a random 
cell of the system. 


The above assumption is mathematically convenient as it 
eliminates the use of the complicated expression of Eq. dUi in 
computing R. The cost of doing so is that attention is restricted 
only to operational points where each AP is almost surely 
active due to the presence of C-UEs alone. However, this heavy 
loaded system scenario is actually of the most interest for 
introducing a D2D mode to the conventional cellular network 
as a traffic offloading method. In addition, it follows from Eq. 
dUi that P{K > 0) > 0.955 with Xc/Xa > 5, i.e., a very good 
approximation of the requirement of Assumption 2 is achieved 
by a moderate/reasonable value of C-UEs traffic load. 

Remark: In general, the requirement P(iC > 0 ) « 1 
of Assumption 2 can be achieved by an appropriately large 
Ac + (1 — p)Xd, i.e., the constraint on Ac can be relaxed by 
taking into account the cellular load due to D-UEs. However, 
this approach imposes a constraint on the mode selection 
probability p that depends on Ac and Xd, complicating the 
closed form solution of the design problem. 

The following proposition provides i? in a simple closed 
form expression. 

Proposition 2: The average normalized UE rate under Con¬ 
straint 1 and Assumption 2 equals 

R - i?„oD2D (^1 + g)) , (16) 

with i?noD 2 D as in Eq. (fTTI) with P(iT > 0 ) = 1 , and 

^ J — p for overlay D 2 D, 

f[P, q) - I _^^^^_(c,qp^+c3gp^-i) _ ^ underlay D2D, 

(17) 

with 7 = 1 , 2 for probabilistic and distance-based 

mode selection, respectively, ci = ^(1 -f p{9o)), C 2 = 

(l/2)Ad7rr|^,^,^K0o/“, ci = ci(l -f X^/Xd), ca = C 2 (l - 
K 6 »Q^“/(l-fp( 6 'o))), andc3 = C2K9l^°‘{l+Xc/Xd)/{l+p{9Q)). 

Proof: Result follows after some lengthy but tedious 
algebra by substituting Eqs. (|7]i, (|9]l, (fTSl i. and (fTJt with 
P(iT > 0) = 1 into Eq. (fTOl i. ■ 

Proposition 2 introduces the function / which plays an 
important role in system design as it is only via / that the 
D2D mode parameters affect R, i.e., the optimization of R 
w.r.t. (p, q) is equivalent to the optimization of / w.r.t. (p, q). 
In addition, / alone indicates whether the incorporation of 
a D2D mode results in i? > i?noD 2 D or i? < i?„oD 2 D for a 
given Pop, depending on whether it is positive or non-positive, 
respectively. Interestingly, the effect of Pop on / is compactly 
described by the set of the coupled parameters ci, Ca G 
for the overlay, and ci, ca, C 3 G for the underlay D2D 
case. Note that / is independent of Ac for the overlay D2D 
case. 

The above observations significantly simplify the represen¬ 
tation (and investigation) of 7 ?.d 2 d which can be equivalently 
written for the overlay case as 

Rdid 

= {(ci,C 2 ) G : 3 {p,q) G (0,1]^ for which f{p,q) > 0} 

= {(ci,C 2 ) G : max f{p,q)>0}, (18) 
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where the second equality follows by noting that f{p,q) is 
uniformly continuous on (0,1]^. The operational region for 
the underlay D2D case can be similarly described in terms of 

(Cl,C2,C3) G 


instead of employing probabilistic approaches for mode selec¬ 
tion and/or channel access procedures. In addition. Scheme 
3-d achieves its best performance with an underlay D2D 
deployment suggesting that the latter exploits the system 
resources more efficiently than an overlay deployment. 


B. Ordering of D2D Deployment Schemes for Heavy Loaded 
System 

Before proceeding with a detailed investigation of the 
operational regions and optimal mode parameters for each 
of the four D2D deployment schemes considered in Sec. 
II.B, significant insights on their merits can already be made 
at this stage of the analysis. The first is provided by the 
following proposition which shows that the flexibility offered 
by employing both a mode selection and a channel access 
procedure provides much less gains than what would probably 
be expected. 

Proposition 3: There do not exist operational points in 
^D 2 D for which the maximization of R is achieved with 
(p*i Q*) G (0,1)^ when overlay D2D (either with probabilistic 
or distance-based mode selection) or underlay D2D with 
distance-based mode selection are employed. 

Proof: A necessary condition for {p*,q*) G (0,1)^ 


df{p*,q) I 


df{p,q* 


= 0. For the Schemes 


dq l'?=9* — dp IP=P* 

referred to in the proposition, this system of equations is either 
inconsistent or gives a solution for which f{p*,q*) = 0 . ■ 

The above result essentially states that Schemes 4-d and 4- 
p for the overlay case and Scheme 4-d for the underlay case 
need not be considered explicitly as they achieve their optimal 
performance when operating as their simplified versions which 
do not allow for (p, q) G (0,1)^ by definition. For the Schemes 
remaining into consideration, an ordering can be made based 
on the following binary relation. 

Definition 4: Relation X A Y (“X greater than Y”), where 
X and Y belong to the set of D2D deployment schemes 
described in Sec. II.B, implies that the maximum rate provided 
by Scheme X is equal to or greater than the maximum rate 
provided by Scheme Y, for any operational point. 

Note from Eq. (fTsT l that X A Y also implies TI% 2 d — 
^D 2 D’ where Ti^ 2 D’ ^d 2 d D2D operational regions 

of Schemes X and Y, respectively. 

Proposition 4: The set of D2D deployment schemes forms 
a totally ordered set with 

Scheme 3-d ^ Scheme 2 A Scheme 3-p A Scheme 1, 
for the overlay case, and a partially ordered set with 
Scheme 3-d A Scheme 4-p A Scheme 2 A Scheme 1, 


and 

Scheme 3-d ^ Scheme 4-p ^ Scheme 3-p ^ Scheme 1, 

for the underlay case. Scheme 3-d of the underlay D2D case 
is greater than Scheme 3-d of the overlay D2D case. 

Proof: See Appendix C. ■ 

Scheme 3-d is therefore the optimal choice for both overlay 
and underlay D2D deployments, showing the significance of 
mode selection based on instantaneous per link information 


VI. Optimal D2D Mode Parameters and 
Operational Region for Heavy Loaded System 

This section explicitly investigates the optimal D2D mode 
parameters and operational region of each scheme which, 
in addition to being of interest in their own right for the 
system designer/operator, will allow to obtain insights on the 
dependance of system performance on Ad and rd^max, the most 
critical operational parameters related to D2D use/business 
cases. 

A. Scheme 1 (Baseline) 

The baseline scheme is the simplest approach for incorpo¬ 
rating D2D communications and has been routinely employed 
in D2D studies, which provides the incentive to investigate it 
even though it is a special case of the other schemes. The next 
lemma follows directly from Eq. (fTTI i by setting p = q = \. 

Lemma 5: The operational region of Scheme 1 equal^ 

Kid = {ci > e^^j, 
for overlay D2D, and 

Kid = {ci > 

for underlay D2D. 

Lemma 5 provides a simple description of the operational 
region of Scheme 1, allowing the system operator to check 
whether the introduction of a baseline D2D mode can provide 
benefits for a system operational point of interest (the latter 
may be determined by potential D2D-related business/use 
cases). To aid the visual comparison between overlay and 
underlay D2D options and examine the effects of Ad and 
Td.max on performance, Eig. 1 depicts 72 .q2d in terms of 
quantities Ad/Aa and AdIE(7rr^) = (l/2)Ad7rr^ for the test 
case of a = 4 and Oq = —6 dB, the latter value corresponding 
to the minimum SIR supported by current wireless standards 
||2^ . Note that Ad/Aa equals the average number of D-UEs 
within a random cell 0, whereas AdE(7rr^) can be interpreted 
as the average number of D-UE sources that are closer to a 
randomly selected D-UE than its own source. Eor the underlay 
case, the region is parametrized by the average number of 
C-UEs within a random cell, Xc/Xa (overlay performance is 
independent of Ac). 

The expressions of Lemma 5 verify the intuition that 
baseline D2D communications cannot provide benefits with 
excessively large rd, max, although the maximum supported 
I'd. max increases with Ad/Aa. Interestingly, the operational 

^The notation of operational region will discriminate among schemes by use 
of a superscript, with no differentiation between overlay and underlay D2D 
cases as this will be clear from context. For simplicity, the sets describing the 
operational region for overlay D2D will be compactly denoted only in terms 
of constraints on (ci, C 2 ) (as defined in Proposition 2) with the understanding 
that (ci,C 2 ) G and similaiiy for the underlay case. 
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Fig. 1. T^o 2 D overlay and underlay D2D (a = 4, dg = —6 dB). Lines 
depict the boundary and arrows point to the interior of the region, respectively. 



Fig. 2. overlay and underlay D2D (a = 4, Oq = —6 dB). Lines 

depict the boundary and arrows point to the interior of the region, respectively. 


region for overlay D2D case requires a minimum Xd/Xa, even 
when rd^max — >■ 0 (c2 —?> 0 ), as it must always hold ci > 1 => 
Xd/Xa > 1/(1 + p{0o))- This might be surprising as with an 
arbitrarily small D2D link distance all D2D transmissions are 
successful and one might expect that D2D communications 
would therefore be beneficial irrespective of Xd in that case. 
However, for sufficiently small Xd, the portion of bandwidth 
dedicated to overlay D2D communications according to Eq. 
(foi l is so small that the average rate achieved by the reliable 
D2D links is in fact smaller than the one provided by the 
less reliable and time shared, but of much higher bandwidth, 
cellular links. 

In contrast, noting that ci > Xc/Xa by definition and 
Xc/Xa > 1 under Assumption 2, underlay D2D is able to 
accommodate any Xd/Xa for arbitrarily small max (C 2 +C 3 —>■ 
0). However, underlay operation imposes stricter constraints 
on the maximum supported max than overlay operation 
when moderate to large Xd/Xa are considered. This difference 
becomes more pronounced for increasing Xc/Xa with underlay 
^D 2 D ^ asymptotically. Considering D2D use/business 
cases that require operation with large D-UE density and 
D2D link distances, it follows that the overlay option is 
more appropriate when Scheme 1 is considered for D2D 
deployment. 

B. Scheme 2 (D-UEs Employ D2D Mode by Default) 

It is straightforward to show the following lemma by setting 
p = 1 in Eq. (fTTI l. 

Lemma 6: The operational region of Scheme 2 equals 
^D 2 D ~ ^D 2 D ^D 2 D’ 7 ^d 2 d i*' Lemma 5, 

“^020 = {Cl > eC2,C2 > 1} 
for overlay D2D, and 

"^020 = {ci > e(c2 + C3), C2 + C3 > 1 } 
for underlay D2D. The rate maximizing access probability 


equals 

^ < 1 for overlay D2D, 

< 1 for underlay D2D, 

when the system operates in 7 ^d 2 D’ < 7 * = 1 otherwise. 

As expected, Til, 2 D ^d 2 D’ he., introduction of a channel 
access mechanism to the baseline scheme can only result in an 
increase of the operational region. This region augmentation 
is due to the inclusion of 7 ^d 2 d which is depicted in Eig. 
2. As can be seen, Tt^ 2 D include operational points 

with small r^ max (for finite Xc/Xa), as in that case the SIR 
of a D2D link is large enough and use of a channel access 
mechanism with q < 1 incurs performance loss. It can also 
be shown that 'R-q 2 d ^d 2 d ^ ^ (compare Eigs. 1 and 
2), i.e.. Scheme 2 not only expands 7 ?.q 2 o but also improves 
performance in the subset of 7 ?.Q 2 r, corresponding to larger 
rd,max values. 

Regarding the overlay/underlay comparison for Scheme 2, 
it can be easily verified that the maximum R provided by 
the underlay version is greater than the corresponding one 
provided by the overlay version for any operational point, 
which also implies a greater operational region for underlay 
Scheme 2 as can also be seen in Fig. 2. 

C. Scheme 3-p (Probabilistic Mode Selection with D2D Links 
Always Active) 

Proposition 5: The operational region of Scheme 3-p 
equals = ^d 2 d ^dId’ with 'R-\) 2 d Lemma 5, 

*^D 2 D “ ^^ 2 ^^ < Cl < ^ 1 < C 2 < 2 ^ 

U {c2e < Cl, C 2 > 2} 
for overlay D2D, and 

'• -X r - gC2+C3 

^D 2 D = < Cl < ^ _ , C2 < U < Ci,C2 > 1} 
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Fig. 3. overlay and underlay D2D (a = 4, 0o = —6 dB). Lines 

depict the boundary and arrows point to the interior of the region, respectively. 


for underlay D2D. The rate maximizing D2D mode selection 
probability equals 

f ^ < 1 for overlay D2D, 

P = < “5 19 

I < 1 for underlay D2D, 

when the system operates in 7 ^q 2 D’ where x* is the unique 
solution of a:e ““(2 — x) = C 2 IC 1 , 1 < x < min{ 2 ,C 2 }, and 
e“®(l — x) = 0 < a: < min{l,C 2}5 for overlay and 

underlay D2D, respectively. For all other operational points of 
'^D2D’ P — 

Proof: See Appendix D. ■ 

Figure 3 depicts f®'^ overlay and underlay D2D. 

Similar observations as for Scheme 1 also hold here, i.e., 
overlay and underlay regions of Scheme 3-p are partially 
overlapping, with the overlay version supporting larger max 
for moderate to large Xd/Xa, whereas underlay version is 
severely limited with increasing Xc/Xa- 

An interesting observation that directly follows from 
Lemma 6 and Proposition 5 is that, for the overlay case only, 
^D 2 D ~ ^D 2 D’ overlay versions of Schemes 2 and 3-p 
provide the same operational region even though under a com¬ 
pletely different approach. However, as shown in Proposition 
4, Scheme 2 is a better choice in terms of maximum R. For the 
underlay case, a clear-cut relation between Schemes 2 and 3-p 
does not exist, as there are operational points included in 7 ^p 2 D 
and not in 7 ^d 2 d versa. However, it can be shown 

that for the operational points belonging to the operational 
regions of both underlay Scheme 2 and Scheme 3-p, Scheme 
2 provides the largest maximum R. 


D. Scheme 3-d (Distance-Based Mode Selection with D2D 
Links Always Active) 

The results for this scheme are presented separately for its 
overlay and underlay versions as the expression of / for the 
latter does not allow for exact closed form expressions (as was 
the case for the schemes examined so far) and necessitates the 
use of bounding techniques. 


Proposition 6: The operational region of Scheme 3-d for 
overlay D2D equals 7 ^o 2 d ~ ^d 2 d ^d 2 D’ with ttt 

Lemma 5 and 


7^, 


3-d 


= < sfficffe < Cl < 


'D 2 D ^ V '-'.L 0/^1 \ 

2(1-C 2 ) 

U {ci > \/2c2e, C2 > 1}. 


1 

2 <C 2 < 


( 20 ) 


The optimal D2D mode selection probability equals 



when operating in where x* is the unique solution of 

ffxe~^{\ — x) = y/oili^ci), 1/2 < a; < min{l,C2}, and 
p* = 1 otherwise. 

Proof: Proof is similar to the case of Scheme 3-p and is 
omitted. ■ 

Proposition 7: The operational region of Scheme 3-d for 
underlay D2D equals 




3-d 

D2D 


= {Cl > 1 }. 


The region for which p* < 1 is bounded as 

^D 2 D — ^D 2 D ’ where 


3-d,(1) ^ 
■D2D — 




3-d. (/3) 
D2D 


= U < Cl < 


gC2-|-C3 

1-/3(c2+C3)’^^"^'^^ ^ /3 


U / Cl > 1,C2 + C3 > 


1 


( 22 ) 


for /3 = 1 , 2 . 

Proof: See Appendix E. ■ 

Proposition 7 provides ttt terms of a simple constraint, 
namely, ci > 1, which, as was stated in the analysis of Scheme 
1, holds for any operational point under Assumption 2, leading 
to the following conclusion. 

Corollary 1: Underlay Scheme 3-d enhances cellular sys¬ 
tem performance in terms of R for any operational point. 

The above result is of signihcant importance as it shows that 
exploitation of knowledge of link distances between D-UEs and 
their sources with an underlay D2D deployment is sufficient to 
achieve optimal performance w.r.t. D2D operational region. In 
this respect, incorporation of other/additional information, e.g., 
distance of D-UEs from their closest AP ifTSll . is not necessary 
(although it may lead to a greater maximum R). 

Eigure 4 depicts overlay D2D and the correspond¬ 

ing region bounds for the underlay D2D case. Note that the 
bounds are tight, with the operational points not included in 
'^D 2 D corresponding to small r^.max values. Eor the overlay 
case, exploitation of D2D link distance information provides 
signihcant region enlargement (compare with previous hgures) 
although there exist operational points corresponding to small 
Xd/Xa and large rd^max that are not included. 


E. Scheme 4-p (Probabilistic Mode Selection and Channel 
Access for D2D Links) 

According to Proposition 3, Scheme 4-p for underlay D2D, 
is the only candidate scheme for which there exist points 
within the D2D operational region achieving maximum R with 
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Fig. 4. ^d 2 D probabilistic mode selection for overlay and underlay 
D2D (a = A, Oq = —6 dB). Lines depict the boundary and arrows point to 
the interior of the region, respectively. 


(P) q) S (0,1)^. This is indeed the case, as described by the 
following proposition whose proof is omitted. 

Proposition 8: The operational region of Scheme 4-p for 
underlay D2D equals 7 ^q2d ~ ^d 2 d ^d 2 D’ 

^D 2 D’ ^D 2 D Lemma 6 and Proposition 5, respectively, 
and 


7^q2d = max 


— ,ec3 < Cl < 
C3 


e(C2 + C3)^ 


C3 


, C2 + C3 > 11 . 

(23) 


The optimal D2D mode parameters equal 


P 


* 


q 


* 



< 1 , 


(24) 

(25) 


when the system operates in p* = 1 and q* as in Lemma 

6 when the system operates in 7 ^q2d \ ^d 2 D’ 9* = 1 ^nd 

p* as in Proposition 5, otherwise. 

The region corresponding to operational points for 

which {p*,q*) G (0,1)^ is shown in Fig. 5. As can be 
seen, the region expansion offered by Scheme 4-p is rather 
limited, with —>■ 0 for increasing Xc/Xa- Note that, 

even though underlay Scheme 4-p allows for both non-trivial 
mode selection and D2D channel access procedures, its region 
does not include all possible operational points which can be 
attributed to the purely probabilistic nature of both schemes 
that do not exploit per link information as in the case of 
underlay Scheme 3-d. 

Having examined all schemes. Fig. 6 compares the oper¬ 
ational regions of overlay and underlay D2D deployments 
achieved with either probabilistic or distance-based mode 
selection. In addition to the already established fact that 
underlay D2D with Scheme 3-d is the best option that can 
support any operational point, it can be seen that an overlay 
deployment with probabilistic mode selection provides the 
smaller operational region, whereas the regions of overlay 



0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 


(l/2)Arf7rr|,„„ 


Fig. 5. for underlay D2D (a = 4, Sq = ~6 dB). Lines depict the 

boundary and~ an'ows point to the interior of the region, respectively. 



(l/2)Aarrr^,„ax 


Fig. 6. Operational regions of overlay and underlay D2D (a = A, 9 q = —6 
dB, Xc/^a = 10). Thick Lines depict the boundary and arrows point to the 
interior of the region, respectively. Thin lines represent the maximum-gain 
level sets of overlay D2D deployment with probabilistic mode selection. 


D2D with distance-based mode selection and underlay D2D 
with probabilistic mode selection partially overlap. With small 
operational Vd^max, probabilistic mode selection may be pre¬ 
ferred as it provides benefits without the cost associated with 
obtaining distance information, however, distance-based mode 
selection always provides larger rates (Proposition 5) and is 
mandatory for large operational max- 

Also shown as thin dashed lines in Fig. 6 for the case of 
overlay D2D deployment with probabilistic mode selection, 
are the maximum-gain level sets, i.e., operational points where 
the ratio i?/i?noD 2 D under optimal (p, q) is constant (each level 
set is labeled by the value of achieved gain). The level curves 
of the other schemes show similar trends and are omitted for 
clarity. It can be seen that incorporation of D2D communica¬ 
tions in the system provides largest gains for operational points 
corresponding to small D2D links and large D-UE density, as 
expected. Note that small perturbations (e.g., due to changes 
in D-UE density) of operational points corresponding to large 
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(a) 


(b) 


Fig. 7. Dependance on normalized ^ax of (a) rate improvement of D2D-enabled system and (b) normalized threshold distance employed by Scheme 3-d. 
Lines depict analytical results and markers depict Monte Carlo evaluations (a = 4, dg = —6 dB, Ac/A^ = A^/Aa = 10). 


D2D link distances result in small changes of gain, whereas 
gain is highly sensitive with perturbations of operational points 
corresponding to small D2D link distances. 

VII. Performance Dependance on Maximum D2D 
Link Distance ans D-UEs Density 

This section further investigates the effect on system per¬ 
formance of Td max and Xd, this time w.r.t. the performance 
gain provided by the introduction of a D2D mode in terms 
of maximum R, denoted as R*. In addition, the presented 
test cases will also serve to validate the accuracy of the 
analytical expressions of the previous sections by comparison 
with Monte Carlo evaluation of Eq. (fTOl) . The latter was 
obtained under the system model of Sec. II under Assumption 
1, i.e., the cellular uplink was not considered. In addition, the 
optimal D2D mode parameters of each scheme were employed 
for each operational point as given in Sec. VI {p* was obtained 
by numerical search for the case of underlay Scheme 3-d). 
In all cases, test values of a = 4 and Oq = —6 dB were 
considered, with Monte Carlo results obtained by averaging 
over 10® independent system realizations, each with 30 APs on 
average. Performance of underlay Scheme 4-p is not depicted 
as it provides insignificant gains compared to Scheme 2 for 
the test cases considered. 

1) Dependance on rd,max' Eigure 7a shows the gain ratio 
^*/^noD 2 D provided by the incorporation of a D2D mode 
to the conventional cellular network as a function of the 
normalized maximum D2D link distance Vd^max'^V^^ where 
l/(2v%r) is the average distance of a UE from its closest 
AP Q. This normalization is convenient, not only because it 
relates the distances involved in setting up cellular and D2D 
links, but also since the achieved rate R depends only on the 
ratios of densities, Ac /Xa and Xd/Xa, and not on their absolute 
values, as can be directly verified from Eqs. ( fTSl l and ([TTT i. Eor 
the results of Pig. 7a, a test case of X^/Xa = Xd/Xa = 10 was 
considered, i.e., there are on average 10 C-UEs and 10 D- 
UEs within each cell, corresponding to a (future) operational 


scenario where D2D use cases constitute a significant part of 
the system load. 

It can be seen that the analytical results for the overlay case 
match almost exactly the Monte Carlo evaluations, whereas for 
the underlay case, they slightly underestimate performance. 
Similar correspondence of analytical and Monte Carlo evalu¬ 
ation results was observed for all operational points that were 
tested (not shown). 

As expected, the performance gain of all schemes (either 
overlay or underlay) diminishes with increasing rd,max- Eor 
I’d,max — 0, all schemes boil down to Scheme 1, with 
R*/RnoD 2 D ~ 12.82,6.67, for underlay and overlay D2D, 
respectively (not shown in the figure). This superiority of 
underlay D2D holds for increasing rd,max values up to 
rd,max ~ 0.6 X l/(2v^). Eor larger rd,max, the ordering of 
the overlay/underlay schemes becomes involved and depen¬ 
dent on rd, max- 

In accordance with the analysis of the previous section, all 
schemes with the exception of Scheme 3-d provide gains up 
to a certain maximum rd,max- An important observation is 
that these maximum values are of the order of the distance 
from the closest AP, i.e., introduction of D2D communications 
can enhance cellular network performance even for distances 
that are not restricted to the common notion of proximal 
communications assuming distance significantly smaller that 
the distance from the cellular APs. 

Regarding Scheme 3-d, it should be emphasized that even 
though it is able to provide gains for large values of rd,max, 
it does so by only allowing the establishment of D2D links 
of distance th = ■\/P'rd, max, which is strictly smaller 

than I’d,max for p < 1. Eigure 7b shows the (normalized) 
threshold value rd,* when p — p* as a function of the 
(normalized) maximum D2D link distance rd.max (curves 
labeled “maximum R”). In particular, for rd,max equal to 
1.5 times the average distance from the closest AP, underlay 
Scheme 3-d only permits establishment of D2D links of 
distances up to about 0.45 times the average distance from the 
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closest AP, with the overlay version being less conservative. 
Note that if larger values of r^^th are needed, e.g., due to 
application requirements, these can be obtained by allowing 
some performance degradation w.r.t. R. In particular, r^^th may 
be set to any positive value not exceeding sup(p)rd_max with 
sup(p) the supremum of p for which R > i?noD 2 D holds. The 
latter can be found numerically and the corresponding bound 
on rd,th is shown in Fig. 7b (curves labeled ‘Vd,* bound”). It 
can be seen that relaxing the performance gain provided by 
D2D transmissions can significantly enlarge th, however, 
still remaining strictly smaller than max when the latter 
exceeds a certain threshold. 

2) Dependance on Xd/^a- Figure 8 shows i?*/i?noD 2 D 
as a function of Xd/Xa for Xc/Xa = 10 , and r^^ max = 
0.8 X l/(2i/A^). This study corresponds to a scenario of 
interest for network operators, where the network is called 
to accommodate for an increasing number of D-UEs due to 
the increased penetration of D2D use cases and introduction 
of D2D enabled devices. 

Note that for this scenario, performance of both con¬ 
ventional and D2D-enabled cellular network decreases with 
increasing Xd/Xa as there are more UEs competing for the 
same set of resources. However, as can be seen from Eig. 
8, i?*/i?noD2D is an increasing function of Xd/Xa when the 
latter is small, i.e., D2D-enabled cellular network performance 
degrades slower than the conventional network since it exploits 
the direct communication possibilities. Eor the overlay case, 
this advantage of D2D-enabled network is observed only above 
a threshold on Xd/Xa (see corresponding discussion following 
Lemma 5). Interestingly, Schemes I and 3-p provide increasing 
performance up to a certain value of Xd/Xa above which 
the performance degrades due to the increasingly high inter¬ 
ference from the D2D transmissions that cannot be handled 
by these two schemes. In contrast, both schemes 2 and 3-d 
show a monotonic ally increasing performance gain as Xd/Xa 
increases. This trend can be verified analytically by direct 
substitution of the expression for q* as provided in Lemma 
5 into in Eq. (m and the fact that Scheme 3-d ^ Scheme 
2. Intuitively, Scheme 3-d provides the largest gains since the 
increased number of D-UEs in conjunction with the distance- 
based mode selection procedure results in establishing many 
D2D links of small link distances, essentially exploiting a form 
of multiuser diversity. 

VIII. Conclusion 

This paper investigated and characterized the operational 
region of D2D communications for enhancing the performance 
of the conventional cellular network in terms of average 
user rate. Various D2D deployment schemes were examined 
for both overlay and underlay options. Eor the important 
case of a heavy loaded network, the optimal D2D mode 
parameters as well as the operational region of every scheme 
were obtained analytically. It was shown that exploitation of 
D2D link distance for mode selection is alone sufficient to 
render D2D communications beneficial for any operational 
point of interest. Under the appropriate D2D deployment 
scheme, significant performance gains can be achieved even 



Fig. 8. Performance gain of the various D2D deployment schemes as a 
function of A^j / Aa . Lines depict analytical results and markers depict Monte 
Carlo evaluations (a = 4 , 0 o = —6 dB, Ac/Aa = 10 ). 

with significant number of D2D-enabled devices and D2D 
link distances of the order of the distance of UEs to their 
closest AP. These observations suggest that the introduction 
of a D2D mode as an add-on to the conventional cellular 
network has significant potential for success as the operational 
points where performance enhancement is observed are not 
restricted to the typically-considered regime corresponding to 
small distances compared to the distance from the closest AP 
and user densities. 

Appendix 

A. Proof of Lemma 1 

Consider the underlay D2D case as the overlay case fol¬ 
lows without taking into account the effect of D2D-generated 
interference. Derivation follows the same steps as in iini 
and is briefly described here as there are a few details that 
are different due to system model differences. Let /ceii = 

SyeiAyo P<^9v\y\~'"^ ^d 2 d = \y “ Xy\°‘gy\y\-‘^ with 

all quantities as defined in Eq. (uJ. It follows from Eq. ([T]) 
and the independence of Idid, that P(SIR > 0) = 
^[J^icc,Myo\“/Pa)Ci^^^{0\yo\°'/Pa)], where £ 3 ,( 5 ) denotes 
the Laplace transform of variable x and the expectation is 
over |j/o|, which is Rayleigh distributed with mean l/{/2.y/X//) 
Q. Approximating as generated from a thinning of 
with retention probability P(Ar > 0) for each x G ^a iflTl . it 
follows that Cl 

{0\yor/Pa) - exp {-7tF{K > O)A,p(0)|yoH • 

It is easy to see, based on the system model and fundamental 
properties of HPPPs |I5], ifT^ . that ^d forms a HPPP of 
density qpXd with each node transmitting with an independent, 
identically distributed (i.i.d.) power Pd, equal to and 
’’'d^i'f'd < Pd.th) under probabilistic and distance-based mode 
selection, respectively. Starting from the Laplace transform of 
the interference power generated by an HPPP with i.i.d. node 
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powers GSl, 



= exp 


= exp 


^-gpAdK7rE(P^^“) 



with 7 = 1, 2 for probabilistic and distance-based mode 
selection, respectively, where the second equality follows by 
noting that E(r^) = (l/2)r^ and E(rp(rd < rd.*)) = 
(l/2)r^jjj = (l/2)pr^ The result of (|2]l then follows by 
evaluating the expectation over ji/ol- 


B. Proof of Lemma 3 

Normalize w.l.o.g. the unit area so that the density of APs 
becomes = 1. By the properties of HPPPs Q, the positions 
of UEs employing cellular transmissions constitute an HPPP 
of density A^ = (Ac + (1 - p)Xd)/K- Let fviy), fxix) 
denote the probability density functions of the area Y of 
the cell containing the origin and the area X of any other 
randomly selected cell, respectively. It can be shown that 
friy) = yfxiy) lUll- Then, 


E 


Ko + 1 


E 




—/(K„ = k) 


OO 

(a) 1 
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(b) J_ 
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iKy)\-x 


c JO 


\k^0 


e ^'^dy 

{Ky)Y.\ 

{k + i)\ 


dy 
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1 - 


fx{y) 


.-Kv 


dy, 


where (a) is due to Kq being Poisson distributed with mean 
X'^y given Y = y, and (b) follows by interchanging integral 
and summation (Fubini’s theorem). Noting that the last integral 
equals P(Ar = 0) completes the proof. 


C. Proof of Proposition 4 

Consider the overlay D2D case. Starting from / provided 
by Scheme 1, the following inequalities hold for any opera¬ 
tional point (ci,C2) G (not necessarily included in the 
operational region of a scheme). 

/(l,l) = cie-^^-l 

< max — p} (26) 

pe(o,i] 

= max {p(cipe“°^P — 1)} 
pe(o,i] 

< max — 1} (27) 

pe(o,i] 

= max — 1} 

pe(o.i] 

< max —p}- (28) 

pe(o,i] 


Noting that the right hand side of Eqs. (l26l l. dZTl l. (l28T l 
correspond to the maximum / provided by Schemes 3-p, 2, 
and 3-d, respectively, completes the proof. The orderings for 
the underlay case and the superiority of the underlay version 
of Scheme 3-d over its overlay version are shown similarly. 


D. Proof of Proposition 5 

Consider the overlay D2D case. For the operational points 
where p* < 1, it must hold = 0, which is 

equivalent to 

xe-^2-x) = ^, (29) 

cl 

where x = p*C 2 < 02 - Since the right hand side of ( l29b is pos¬ 
itive it must hold x < min{2,C2}. Assuming this is the case, 
f{p*,0 = f{x/c 2 ,l) = which requires x > 1 in 

order to have a positive value, resulting in x S (1, min{2, C2}), 
which in turn requires C2 > 1. Assuming the latter holds, and 
noting that xe““(2 — x) is a strictly decreasing function of 
X in 1 < X < 2, it is easy to see that a unique solution 
of Eq. (|29] | exists as long as C2/C1 > 1/e, for C2 > 2 and 
1/e < C2/C1 < C2e“'^^(2 — C2) for 1 < C2 < 2, which 
provide the description for 7 ^d 2D' Derivation of the result for 
the underlay case is similar. 


E. Proof of Proposition 7 

Clearly, 7 ^d 2 d = ^d 2 d D ’ 7 ^d 2D’ where 7 ^d 2 d ’^^e set 

of operational points for which p* < 1 and 71^20 
Lemma 5. Noting that the maximum / for underlay Scheme 
3-d can be bounded as maXpg(Q — p} < 

maXpg(o,i)/(p, 1) < maxpg(07){cipe“(‘==+^3)p - p}^ it 

follows that ^ ^D2D C with 

/3 = 1, 2, denoting the sets of operational points where 
the lower and upper bounds are positive, respectively. The 
bounding sets can be obtained as in Eq. (l22l l using a similar 
procedure as previous derivations and it is easy to see that 
Kxd U 7 ^d 2 d''^ = {ci > 1} for ^ = 1, 2. 
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